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ABSTRACT 

We present an analysis of the mid-infrared emission lines for a sample of 12 low metallicity Blue Compact 
Dwarf (BCD) galaxies based on high resolution observations obtained with Infrared Spectrograph on board the 
Spitzer Space Telescope. We compare our sample with a local sample of typical starburst galaxies and active 
galactic nuclei (AGNs), to study the ionization field of starbursts over a broad range of physical parameters 
and examine its difference from the one produced by AGN. The high-ionization line [OIV]25.89/^m is detected 
in most of the BCDs, starbursts, and AGNs in our sample. We propose a diagnostic diagram of the line 
ratios [OIV]25.89^m/[SIII]33.48/im as a function of [NeIII]15.56Aim/[NeII]12.8lAim which can be useful in 
identifying the principal excitation mechanism in a galaxy. Galaxies in this diagram split naturally into two 
branches. Classic AGNs as well as starburst galaxies with an AGN component populate the upper branch, 
with stronger AGNs displaying higher [NeIII]/[NeII] ratios. BCDs and pure starbursts are located in the lower 
branch. We find that overall the placement of galaxies on this diagram correlates well with their corresponding 
locations in the log([NII]/Ha) vs. log([OIII]/H/?) diagnostic diagram, which has been widely used in the 
optical. The two diagrams provide consistent classifications of the excitation mechanism in a galaxy. On the 
other hand, the diagram of [Nelll] 1 5.56/im/[NeII] 12.8 l^rnvs. [SIV] 10.5 1 Aim/[SIII] 1 8.7 1 /^m is not as efficient 
in separating AGNs from BCDs and pure starbursts. Our analysis demonstrates that BCDs in general do display 
higher [NeIII]/[NeII] and [SIV]/[SIII] line ratios than starbursts, with some reaching values even higher than 
those found at the centers of AGNs. Despite their hard radiation field though, no [NeV] 14.32/xm emission has 
been detected in the BCDs of our sample. 

Subject headings: dust, extinction - infrared: galaxies - galaxies: dwarf - galaxies: starburst - galaxies: active 
- galaxies: abundances 
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1. INTRODUCTION 

Starburst galaxies are hosts of intense star formation activ- 
ity with star formation rates (SFR) more than a factor of 3 
compai-ed to the Milky Way (SFR>10Moyr"'). The stellar 
population of the starbursts, especially the number of mas- 
sive stars and the slope of the initial mass function (IMF), 
directly affect the hardness of the ionization field in the in- 
terstellar medium (ISM). Therefore, ratios of emission lines 
of different excitation levels are powerful tools to study the 
properties of the ionization field and to infer the stellar pop- 
ulation of the starbursts. Infrared emission-line ratios are of 
particular interest and have been explored extensively with the 
Infrared Space Observatory (i.e. Thornley et al. 2000; Forster 
Schreiber et al. 2001; Verma et al. 2003, 2005), as the active 
star forming regions in the disks of starburst galaxies and cir- 
cumnuclear regions are often optically obscured by dust. As a 
result, the emitted radiation, originating from a massive star- 
burst and/or an accreting active nucleus, is emerging in the 
mid-infrared (mid-IR) (see Vacca et al. 2002; Gorjian et al. 
2001) and far-IR. The advent of the Spitzer Space Telescope 
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with its improved sensitivity opens new avenues in this direc- 
tion, enabling detailed studies of large samples of quiescent 
and star forming galaxies (Dale et al. 2006, 2009; O'Halloran 
et al. 2008; Bernard-Salas et al. 2009), as well as more com- 
plex dust enshrouded infrared luminous sources which often 
host an active galactic nucleus (AGN) (Armus et al. 2007; 
FaiTah et al. 2007). 

Estimating the hardness of the ionization field plays a key 
role in determining the dominant power source at the centers 
of galaxies. Nuclear starbursts and accreting super-massive 
black holes (SMBH) are often intertwined and coexist, as best 
demonstrated in the case of Ultra Luminous Infrared Galax- 
ies (ULIRGs)(see e.g. Genzel et al. 1998; Armus et al. 2007, 
and references therein). The radiation of an AGN, produced 
from the accretion of gas and stars onto the SMBH, has a 
much harder spectrum compared to the one produced by the 
massive young stellar populations. Therefore, the hardness 
of the radiation field has been a key parameter addressed by 
previous diagnostic methods to determine the energy source 
in galaxies (e.g. Baldwin et al. 1981; Veilleux & Osterbrock 
1987; Genzel et al. 1998; Laurent et al. 2000; Nardini et al. 
2008). 

In the study of the ionization field of galaxies with massive 
starbursts, the case of low-metallicity Blue Compact Dwarfs 
(BCDs) is particularly interesting. BCDs are galaxies char- 
acterized by their blue optical colors (Mb <- 18 mag), small 
linear sizes (<lkpc) and low luminosities (Thuan & Martin 
1981; Gil de Paz et al. 2003). Even though BCDs are defined 
based on their morphological properties, they have been ubiq- 
uitously found to have low oxygen abundances as measured 
from the ionized gas in HII regions. 
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It is well established that the hardness of the stellar radia- 
tion increases with decreasing metallicity (e.g. Campbell et al. 
1986). As a result, the ionization field of BCDs is generally 
much harder compared to a typical starburst. It is also known 
that the mass of a super massive black hole (SMBH) is cor- 
related with the mass of the bulge of the host galaxy (Magor- 
rian et al. 1998; Gebhardt et al. 2000; Ferrarese & Merritt 
2000) while only a fraction of the most massive SMBH are 
actively emitting radiation as a result of mass accretion (Hao 
et al. 2005b). BCDs in general are not massive enough to host 
such an active SMBH, with possibly only a few exceptions 
(see Izotov et al. 2007; Izotov & Thuan 2008). Consequently, 
the excitation mechanism for the hard radiation field in BCDs 
is likely significantly different from the one responsible for 
the radiation field produced by an AGN. Furthermore, unhke 
starburst galaxies which are often obscured by dust, an anal- 
ysis of a large sample of BCDs suggests that in BCDs nearly 
always the young massive stars which ionize the interstellar 
medium are easily detected and their dust extinction measured 
from the optical Balmer fines is typicaUy fairly low (Wu et al. 
2008b). There might be exceptions indicated from the mid-IR 
observations, with SBS0335-052E as the most extreme case. 
However, this is rather a unique galaxy and deviates from the 
usual trends seen in BCDs for different reasons as discussed 
in detail by Houck et al. (2004b). Understanding the differ- 
ences in the excitation mechanism between BCDs and AGNs 
can provide insightful clues to decipher the power source in a 
starburst galaxy also hosting an AGN. Identifying those dif- 
ferences from the mid-IR spectral features is one of the main 
motivations behind the present work. 

BCDs are often much fainter than typical starburst galaxies 
in the infrared. Up until the launch of the Spitzer Space Tele- 
scope (Werner et al. 2004), studies of starburst galaxies in the 
IR have included only a handful of BCDs (Lutz et al. 1998; 
Thornley et al. 2000; Verma et al. 2003; Madden et al. 2006). 
The superb sensitivity of the Spitzer allows us to observe a 
large sample of BCDs (see Wu et al. 2006, 2008a). In this pa- 
per, we provide the ionized gas diagnostic analysis of the low 
metallicity BCDs accessible to the high-resolution modules 
of the Infrared Spectrograph (IRS, Houck et al. 2004a). In 
addition, we contrast the BCD sample with typical starbursts 
and AGNs, to probe the excitation mechanism in these types 
of objects and understand where the difference comes from. 

In §2, we introduce the sample and the data reduction 
method. In §3, we study the high-ionization emission fines, 
and propose a new diagnostic diagram to determine the en- 
ergy source in a galaxy. In §4, we focus on other infrared 
emission lines and compare various excitation diagnostic dia- 
grams of starbursts and BCDs. Our findings are summarized 
in §5. 

2. SAMPLE AND DATA REDUCTION 

2.1. The BCD sample 

We have compiled a sample of 12 BCDs and observed them 
in starring mode with the high-resolution modules (short- 
high, SH: 4.7"xll.3" and long-high, LH: ll.l"x22.3") of 
the IRS as part of the instrument team's Guaranteed Time Ob- 
servation (GTO) programs. The galaxies were drawn from the 
sample of Wu et al. (2006). We exclude Toll214-277 as it is 
too faint and its high-resolution spectrum has low signal-to- 
noise ratio. We also exclude CG0752, CG0563, and CG0598 
as they do not show typical BCD signatures. The characteris- 
tics of individual objects in the sample and analysis based on 
their Spitzer low-resolution spectra are presented in Wu et al. 



(2006, 2008a,b). 

The high-resolution module slits were centered on the nu- 
cleus of each galaxy, as was the case of the low-resolution 
observations (Wu et al. 2006). The data were processed by 
the Spitzer Science Center (SSC) data reduction pipeline ver- 
sion 13.2. Individual pointings to each nod position of the 
slit were co-added. The median of the combined images were 
extracted using the full slit extraction method of the IRS data 
analysis package SMART (Higdon et al. 2004). Then the ex- 
tracted spectra were flux cafibrated by multiplying tfiem witfi 
a relative spectral response function (RSRF), which was cre- 
ated from the IRS standard star, ^ Dra, for which accurate 
templates were available (Cohen et al. 2003). Note that emis- 
sion from the sky background has not been subtracted from 
our spectra due to the lack of dedicated sky observations, thus 
discontinuity exists in some source spectra between SH and 
LH, however, this will not affect the results of our analysis 
since our conclusions are based on ratios of emission lines 
observed in the same module. 

In the Appendix, we include close up plots of strong and 
weak emission lines visible in the IRS high-resolution spec- 
tra of BCDs. In order to measure the line strengths, a Gaus- 
sian profile is fitted to the lines above a local continuum fit- 
ted by an one-order polynomial. The [OIV]25.89/Ltm and 
[FeII]25.99/im are close to each other and they are fit to- 
gether using a two-Gaussian profile with fixed rest-frame 
wavelengths. The measured line strengths for the BCD sam- 
ple are listed in Table 1. The errors are estimated based on the 
error propagation of the Gaussian fit assuming the error of the 
flux as the root mean square {rms) of the continuum fit. Note 
that this does not take into account all the systematic uncer- 
tainties, especially the flux calibration. We estimate these to 
be 5% and choose it as our minimum uncertainty. The upper- 
limits are obtained by measuring the flux of a Gaussian with 
a height three times the local rms, and with a FWHM equal to 
the instrumental resolution. 

No extinction correction has been performed in our anal- 
ysis. The corrections are typically minor for mid-IR fine- 
structure line ratios, except for ratios involving lines found 
close to the silicate features, such as the [SIV] 10.5 1/xm line. 
In the analysis of diagrams involving the [SIV] fine (see Fig- 
ure 9), we provide arrows indicating the effects of an Av=30 
mag of extinction, assuming a uniform dust screen and a Li & 
Draine (2001) extinction curve (see Farrah et al. 2007). 

2.2. The Starburst and AGN sample 

In order to study the ionization field of starbursts over a 
broad range of physical properties, we also examine a sam- 
ple of 24 local starbursts observed in staring mode with the 
IRS high-resolution modules. The details of the sample, their 
high-resolution spectra and fine measurements are described 
in Bernard-Salas et al. (2009). It should be emphasized that 
even though we broadly use the term "starburst" when refer- 
ring to them, some may also have a weak AGN component. 

In addition, we created a comparison sample of typical 
AGNs to investigate how the difference in the excitation 
mechanism of the gas, namely the radiation produced by an 
accretion disk versus massive young star formation, is re- 
flected in the mid-IR emission lines. It is not our intention 
to have an AGN sample that is complete, but rather one that is 
representative of typical AGNs. Hence, we choose bona-fide 
AGNs, with no ambiguity in their classification, where the 
nuclear emission is known to dominate the integrated spec- 
tra of their host galaxies. We select from the Type 1 sources 
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(Quasars and Seyfert Is) used in Hao et al. (2007) for which 
IRS high-resolution observations are available. We also ex- 
amine the strength of the Polycyclic Aromatic Hydrocarbon 
(PAH) emission in their mid-IR spectra and exclude sources 
with strong PAH features, since this would suggest a signifi- 
cant contribution from star forming regions in the disk of the 
host galaxy. Among those sources with weak PAH emission, 
we choose 10 AGNs that show strong mid-IR emission lines, 
in order to obtain the best S/N possible in our measurements. 
The AGN IRS spectra were extracted and the emission lines 
measured in the same way as the BCDs (see §2.1). 

Furthermore, whenever applicable, we also compare the 
properties of emission lines in BCDs with those in low- 
redshift ULIRGs as well as HII regions found in our Galaxy, 
and in the Large and Small Magellanic Clouds (LMC/SMC). 
The sample and line measurements of the low-redshift 
ULIRGs are taken from Farrah et al. (2007), and the HII re- 
gion data from Peeters et al. (2002) and Lebouteiller et al. 
(2008). 

3. fflGH-IONEATION EMISSION LINES DIAGNOSTICS 
3.1. The [OIVJ25. 89 nm Emission 

The IRS high-resolution spectra of BCDs are rich in 
fine-structure emission lines. All BCDs display strong 
[NeIII]15.56/im, [SIII]18.71;um, and [SIV]10.51/im emis- 
sion. The [Nell] 12.8 1 fim Une is well detected in about half of 
the BCD sample, but it is not seen in SBS0335-052E, Tol65, 
and VIIZw403. In the comparison samples of starbursts and 
AGNs, all these lines are strong and we discuss their relative 
flux ratios in §4. 

The presence of high-ionization emission lines is of great 
interest, since their excitation often requires extreme condi- 
tions. One particular case is the [OIV]25.89/imline which has 
an ionization potential (IP) of 54eV, just above the Hell edge, 
and therefore it cannot be readily generated in the HII regions 
surrounding main- sequence stars. However, the [OIV] line 
had already been detected in several starburst galaxies with 
ISO (see Lutz et al. 1998). A possible explanation for this 
was that the [OIV] emission is due to a dust enshrouded AGN 
as the hard radiation field produced by the AGN can easily 
provide the energy needed to ionize O^"^. However, Lutz et al. 
(1998) showed that at least in some starburst systems (such 
as the outer regions of M82), the [OIV] cannot be attributed 
to a weak AGN. Instead it must be due to either very hot 
stars, such as the Wolf-Rayet (WR) stars (see also Schaerer 
& Stasifiska 1999) or to the ionizing shocks associated with 
the starburst activity. 

With the superb sensitivity of the Spitzer Space Telescope, 
we now find that the [OIV] line is fairly common in our BCD 
sample as well as the comparison sample of starburst galaxies. 
We detected [OIV] emission in 8 of the 12 BCDs, and in 21 of 
the 24 starbursts (see Bemard-Salas et al. 2009). This is very 
interesting since, as we will discuss in the foUowing sessions, 
line ratios involving [OIV] provide a good diagnostic tool for 
identifying the presence of an AGN (see also Genzel et al. 
1998; Armus et al. 2007). 

3.2. The [OIV]/[SIII] vs. [NeIII]/[NeII] diagram 

In Figure 1, we plot the line ratios of 
[OIV]25.89/im/[SIII]33.48/Ltm as a function of 
[NeIII]15.56^m/[NeII]12.81^m for our BCDs, starburst 
galaxies AGNs and ULIRGs. We choose to use the line 
ratio of [OIV]25.89/im and [SIII]33.48;um because of their 
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Fig. 1.— A plot of the [OIV]25.89Atm/[SIII]33.48^tm as a function of 
[NeIII]15.56/im/[NeII]12.81/im, for BCDs (diamonds), starbursts (squares), 
AGNs (crosses), and ULIRGs (plus signs). A circle around a symbol indi- 
cates that [NeV]14.32/im emission has been detected, while a filled symbol 
implies that this galaxy also has been classified as a Wolf-Rayet (WR) in the 
optical (see §3.4). The error bars of the line ratios in this and all subsequent 
plots are similar to the symbol size. The galaxies appear to be locaHzed on 
two branches. One branch consists mainly of AGNs, or starbursts/ULIRGs 
hosting an AGN, while the other mainly contains BCDs or pure starburst 
galaxies. Two dashed lines are drawn as a guide to identify the two branches. 

large difference in the ionization potential, and the absence 
of aperture effects, as both lines reside in the same IRS 
long-high module (see also Figure 3 in Genzel et al. 1998) 
. We use the [NeIII]/[NeII] Une ratio as the other variable, 
as it is commonly used as an indicator of the strength of the 
EUV radiation field (e.g. Lutz et al. 1998; Thornley et al. 
2000; Verma et al. 2003). Since both [NelllJ and [Nell] 
lines are of the same element, the ratio is independent of the 
gas abundance. We clearly see that the [OIV |/| Sill] ratio 
separates AGNs from BCDs, with AGNs having in general 
a factor of 10 higher [OIV]/[SIII] ratios than BCDs, even 
though they have similar [Nelll]/[Nell] ratios. Starburst 
galaxies display smaller [NeIII]/[NeII] ratios than AGNs 
and BCDs, but their values of the [OIV]/[Slll| ratio vary. 
In particular, it appears that starburst galaxies populate two 
branches with a common origin (the dashed lines in Figure 1), 
one pointing towards the locus of AGNs and the other to the 
BCDs. Among the starburst galaxies, Mrk266, NGC1365, 
NGC2623, NGC4194, NGC4945, NGC520, NGC660, 
NGC253, NGC1614, NGC1097, NGC3079, NGC4676, 
and NGC7252 are located in the upper branch reaching to 
the AGN locus, while NGC1222, NGC2146, NGC3256, 
NGC3310, NGC3628, NGC4088, and NGC7714 are in the 
lower branch that ties to the BCDs. The galaxies IC342, 
NGC3556, and NGC4818 have no detected [OIV] emission 
and their upper limits are indicated with arrows on the plot. 

Given that the gas in AGNs and BCDs is excited by dis- 
tinctly different physical mechanisms, the fact that they are 
found at the corresponding ends of the two branches of the di- 
agram, also suggests that the excitation mechanisms leading 
to the [OIV] line emission in starbursts galaxies populating 
the two branches is different. This could be understood if the 
mid-IR emission of starbursts in the upper branch is due to 
some AGN component, while in starbursts of the lower branch 
is only due to the reprocessing of pure stellar radiation. There- 
fore, such a diagram can be used as a diagnostic tool to deter- 
mine the origin of the dust enshrouded energy source in star 
forming galaxies. Interestingly, among the starbursts located 
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in the upper branch, six also have detected [NeV]14.32/im 
emission (see Bernard-Salas et al. 2009), which has an IP of 
97eV and is a clear AGN indicator. Five galaxies (NGC520, 
NGC253, NGC1614, NGC4676 and NGC7252) which do not 
have a confirmed AGN signature so far, are located in the 
lower left corner of the plot, where the two branches merge. 
Since at this area of the plot classification is more ambiguous, 
it is possible that a weak AGN component could be present, 
even though its [NeV] emission is too weak to be detected. 

Unfortunately, in most of the ULIRGs with available IRS 
high resolution spectroscopy, their [SIIIJ33.48^m emission 
lines are redshifted out of the wavelength coverage of the IRS 
(see Farrah et al. 2007). Among the 10 ULIRGs that have both 
[OIV]25.89/im and [SIII]33.48/im detected, 9 are located on 
the upper branch and they all have detected [NeV] emission. 
IRAS09022-3615 is the only ULIRG with no [NeV | detection 
and it is located between the two branches. This is again con- 
sistent with the assertion that star forming galaxies with an 
AGN component do follow the upper branch. 

The lower branch of the diagram seems to be populated by 
sources in which the ionizing radiation is due only to stel- 
lar photospheres. In addition to the BCDs, the pure starburst 
galaxies NGC7714, NGC3310, and NGC1222 (e.g. Brandl 
et al. 2004; Wehner et al. 2006; Beck et al. 2007) are located 
in the lower branch before the two branches merge. Moving 
along the branch towards higher values of the [NeIII]/[NeII] 
ratio, the metallicity overall decreases (see §4). The harder 
radiation field produced by the low metallicity environment 
strongly increases the [NeIII]/[NeII] ratio (by up to a factor of 
100), but increases only moderately the [OIV]/[SIII] ratio (by 
a factor of 10). This was also noticed by Lutz et al. (1998), 
who indicated that low-metallicity dwarfs exhibit relatively 
stronger [OIVJ than low-excitation starbursts, but the effect 
is not very pronounced compared to the large excitation dif- 
ference as measured by the [NeIII]/[NeII]. One galaxy with a 
detected [NeV] emission (NGC3628, see Bemard-Salas et al. 
2009) is on the lower branch. We note though that its location 
is also at the far left where the two branches merge, and as a 
result it does not contradict our bifurcation scheme. 
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Fig. 2.— A plot of the [OIV]25.89^jm/[SIII]33.48/im line ratio as a func- 
tion of the [SIV] 10.5 l^m/[SIII]l 8.7 l^m for our samples. The symbols are 
the same as in Figure 1 . Note that the bifurcation into two branches, indicated 
again with the dashed lines, is also visible. 



As with the [NeIII]/[NeII] ratio, the ratio of 
[SIV]10.51;tim/[SIII]18.71;tim has also been widely used 
as a probe of the hardness of the ionization field. When 



plotting the [OIV]25.89Aim/[SIII]33.48^m as a function of 
the [SIV]/[SIII] line ratio, we observe a similar bifurcation 
signature (see Figure 2). The sources located on the upper 
and lower branches in Figure 2 are the same to the ones found 
on the upper and lower branches in Figure 1. The separation 
of the two branches though is not as clear as in Figure 1, 
mainly because of the narrower range of the [SIV]/[SIII] 
values for our sample than the range of the [NeIII]/[NeII] 
ratios. In addition, in several starbursts (8 out of the 24 
galaxies) no [SIV] emission was detected, thus making a 
clear definition of the merging points of the two branches 
challenging. 

3.3. [ OIV] Emission as an AGN Diagnostic 
3.3.1. The [OIV]/[SIII] ratio and optical classification 

Close examination of the [OIV]/[SIII] vs. [NeIII]/[NeII] 
diagram presented in the previous section suggests that it 
can be used as an AGN-starburst diagnostic, similar to the 
log([NII]/Ha) vs. Iog([OIII]/H/3) diagram (e.g. the BPT dia- 
gram, Baldwin et al. 1981), which has been used traditionally 
to distinguish AGNs from starbursts in the optical (Veilleux & 
Osterbrock 1987; Kewley et al. 2001; Kauffmannet al. 2003). 
In both diagrams, star forming galaxies are distributed along 
two separate branches, each having one end respectively an- 
chored by the location of AGNs and BCDs. In Figure 3, we 
plot on the BPT diagram ^40,000 galaxies from the Sloan 
Digital Sky Survey (SDSS, see Hao et al. 2005a). Strong 
AGNs have high [NII]/Hq! and |0IIIJ/H/3 ratios and they are 
found at the tip of the right branch. Low-metaUicity BCDs 
have high [OIII]/H/3 but low |NII|/Ha ratios and they are lo- 
cated at the tip of the left branch. The natural separation of 
galaxies into two branches indicates again a different excita- 
tion mechanism. The blue dotted line, empirically defined by 
Kauffmann et al. (2003), separates the two branches. These 
authors classify galaxies below the line as star forming galax- 
ies, and those above it as AGNs. The red dashed line is taken 
from Kewley et al. (2001), and demarcates the maximum po- 
sition that can be obtained by pure photo-ionization models. 
Galaxies located above the line require an additional exci- 
tation mechanism, such as an AGN, or strong shocks. The 
general classification scheme using the two separation lines 
(Kewley et al. 2006) considers objects above Kewley's line as 
AGN dominated sources, between Kewley's and Kauffmann's 
line as composite AGN and starburst sources, and below the 
Kauffmann's line as star forming galaxies. 

To quantitatively check the analogy between our [OIV] di- 
agram (Figure 1) and the log([NII]/Ha) vs. log([OIII]/H/?) 
diagram used in the optical, we have collected from the litera- 
ture existing optical data for our starbursts and BCDs. Special 
care was taken to match the optical apertures to the 11" wide 
Spitzer IRS-LH slit, in which the [OIV] line is located. This 
is especially important for the starburst galaxies as many of 
them are extended. The most appropriate optical datasets for 
the majority of galaxies were the integrated spectra obtained 
by Moustakas & Kennicutt (2006) from which we retrieved 
the optical spectra data for 13 starbursts and 8 BCDs. 

To expand our dataset, we also use the data for star- 
forming galaxies from Dale et al. (2006), which have match- 
able Spitzer high-resolution mapping spectra extracted with 
an aperture of 23"xl5", and optical spectra taken with an 
aperture of 20"x20". There are 19 objects with both de- 
tections in [OIV]25.89/;im and optical lines used in the BPT 
diagram, which are displayed in the Figure 2 of Dale et al. 
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Fig. 3.— The optical log([NII]/Ha) vs. log([OIII]/H/3) diagram in which 
we include ~40,000 SDSS galaxies, marked as black dots. The red dashed 
hne and the blue dotted line are taken from Kewley et al. (2001) and Kauff- 
mann et al. (2003) and they are used to separate AGNs, to the top right of the 
plot, from starburst galaxies to the bottom left. We include on the diagram 
our sample of BCDs and starbursts as well as the star-forming galaxies from 
Dale et al. (2006) (D06). Note that some galaxies in our starburst sample and 
in D06 also host an active nucleus and as a result are found in the AGN locus. 
Galaxies located at the upper and lower branch of Figure 4 are plotted in 
orange and cyan respectively (see text). 
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as stars). As before, filled symbols are galaxies with WR characteristics in the 
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line are blue. The dashed lines are the same as in Figure 1. The dotted line is 
Equation (1), defined empirically as the separatrix between the two branches. 



(2006). In our Figures 3 and 4, we place our samples with 
optical data, as well as the star-forming galaxies from Dale 
et al. (2006), in the corresponding optical log([OIII]/H/3) vs. 
log([NII]/Ha) and infrared [OIV]/[SIII] vs. [NeIII]/[NeII] di- 
agrams. The colors of the data symbols in each diagram are 
assigned by their corresponding classification in the other di- 
agram (see Figure captions). In Figure 4, the galaxies found 
above Kewley's line in Figure 3 are red^, between Kewley's 
and Kauffmann's line are green, and below Kauffmann's line 
are blue. We found that most optically identified strong AGNs 
(marked with red x) and starbursts (blue squares) are well 

' We do not have the matching optical data for the AGNs in our sample, 
but they are expected to lie well to the light of the Kewley's line. 



separated in the upper and lower branches, especially for val- 
ues of the [NeIII]/[NeII] ratio greater than 0.3. There are 
four galaxies optically classified as starbursts (blue squares) 
located in the upper branch in Figure 4. However, we note 
that their location in the optical diagram in Figure 3 is near 
the locus where the two branches merge, an area where clas- 
sification of individual objects is ambiguous. 

In the optical diagram of Figure 3, the data points are col- 
ored based on their classification in the infrared diagram of 
Figure 4. More specifically, we empirically define a line sepa- 
rating the upper and lower branch, having the functional form 
of: 

/og([OIV]/[SIII]) = 0.08 X log([NeIII]/[NeII]) + 0.014 (1) 

Galaxies located on the upper branch above the line are 
plotted in Figure 3 in orange, and below the line in cyan. This 
demonstrates that the upper and lower branch of the infrared 
diagram correspond well to the right and left branch of the 
optical. 

We should stress once more that the optical spectra are col- 
lected from various sources and have been obtained with dif- 
ferent apertures. This likely introduces uncertainties in the 
separatrix mentioned above, which cannot be easily quanti- 
fied. Even for datasets where optical and infrared apertures 
are matched, wavelength dependent extinction will result in 
discrepancies between the different classifications, in partic- 
ular for dusty sources, since the infrared can probe into re- 
gions which are optically thick. Therefore, it is not surpris- 
ing to have some outliers that do not correspond to the same 
regions in both diagrams. For example, galaxies located be- 
tween Kewley's and Kauffmann's line (green objects in Fig- 
ure 4) may lie either on the upper or lower branch of the in- 
frared diagram. Despite these uncertainties, the overall con- 
sistency of the two diagrams is very encouraging, convincing 
us that the [OIV]/[SIII] vs. [NeIII]/[NeII] diagnostic along 
with equation (1) can be used to classify the energy source of 
nuclear emission in galaxies, particularly for those with high 
dust content. 

3.3.2. The [OIV]/[SUI] ratio and other infrared diagnostics 

The ratio of the flux of [OIV]25.89/im line with the flux of 
lower-ionization line has been akeady proposed to diagnose 
the dominant source of power in ULIRGs (Genzel et al. 1998; 
Armus et al. 2007). In the Genzel et al. study, the ratio of 
[OIV]25.89Aimto the [NeII]12.81^m (or [SIII]33.48Aim) was 
plotted against the strength of the 1.1 fim PAH. The PAH fea- 
tures are strong in starbursts and weak in AGNs, as the PAH 
molecules are easily destroyed by the hard radiation of the ac- 
cretion disk or outshined by the hot dust continuum emission 
from the AGN torus. In Figure 5, we plot the [OIV]/[SIII] 
ratio as a function of the equivalent width (EW) of the 6.2^m 
PAH feature. The PAH EWs for the galaxies in our sample 
are collected from Brandl et al. (2006); Wu et al. (2006), and 
Spoon et al. (2007). As was the case with Figure 3 colors of 
the galaxies in Figure 5 are based on their mid-infrared clas- 
sification in Figure 4. Those found on the upper branch above 
the dotted separatrix line are marked in orange, while those 
on the lower branch are in cyan. We find that the two color 
groups define two separate trends, with the AGNs residing at 
the tip of one group and the BCDs at the tip of the other. The 
upper trend is similar to the trend found in the traditional Gen- 
zel et al. diagram (indicated with the upper dashed line), while 
the other trend clearly falls below it. This suggests that the 
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method to diagnose the fraction of AGN and starburst power 
in the mid-IR for a given source suggested by Genzel et al. 
(1998), can only be applied to sources which are part of the 
upper trend in Figure 5. The upper and lower dashed lines can 
be represented as equation (2) and (3) respectively: 

/o^([OIV] / [SIII]) - log(6 .0) = 0. 1 9/(log(EW(6 .2)) - log(0.7)) 

(2) 

Zog([OIV]/[SIII])-log(0.35) = 0.51/(log(EW(6.2))-log(0.7)) 

(3) 
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Fig. 5.— The [OIV]25.89Atm/[SIII]33.48^m ratio as a function of tlie 
6.2fj,m PAH EW. Starburst galaxies (squares), BCDs (diamonds) and AGNs 
(crosses) located at the upper and lower branch of Figure 4 are plotted in 
orange and cyan respectively. Filled symbols are again WR galaxies. Two 
dashed lines are drawn to indicate the trend of the upper and lower branches. 

Dale et al. (2006) have also proposed a diagnostic dia- 
gram based on the [SIII]33.48/xm/[SilI]34.81/im line ratios 
as a function of the [Nelll]15.56/im/[Nell] 12.8 1/im ratios to 
distinguish AGN from starburst emission. Their method re- 
lies on the fact that [Sill], having an IP of only 8.2 eV, is a 
low-ionization emission line and a strong coolant of X-ray- 
irradiated gas (Maloney et al. 1996). AGNs usually have 
stronger [Sill] emission than starburst galaxies and diagnos- 
tic diagrams involving [Sill] present an advantage in cases of 
faint sources where [OIV] is not easily detected. In Figure 6, 
we plot the galaxies from our sample on the [Slll]/[Sill] vs 
[Nelll]/[Nell] diagram, similarly to the Figure 9 of Dale et al. 
(2006). The galaxies are divided into four regions, and are 
classified accordingly in a statistical sense: galaxies located 
in regions 1 and 11 should be classified as an AGN with a Icr 
confidence interval of 83%-97% and 73%-88% respectively, 
and galaxies found in regions 111 and IV should be classified as 
a star-forming systems with a la confidence interval of 84%- 
93% and 91%-98%, respectively. In Figure 6 , the galaxies 
are also colored the same way as in Figure 3 and Figure 5. We 
find that AGNs and star-forming galaxies situated on the up- 
per branch are preferentially located in regions 1 and 11, while 
BCDs and star forming galaxies on the lower branch are found 
in regions 111 and IV. This suggests that the two methods over- 
all agree. 

3.4. The nature of the [OIV]25.89fim emission in starbursts 

As mentioned in the previous section, the [01V]25.89/im 
emission in galaxies cannot originate from the Hll regions 
around main-sequence stars. Instead, it has to be associated 
either with an AGN component, or with WR stars (Crowther 
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Fig. 6.— The diagram of [Sin]33.48/im/[SiII]34.81/im as a function 
of [NeIII]15.56^m/[NeII]12.81/im proposed by Dale et al. (2006). The 
diagram is divided into four regions (I-IV). Galaxies located in regions I 
and II should be classified as AGN with a Icr confidence interval of 83%- 
97% and 73%-88% respectively, and galaxies found in regions III and 
IV should be classified as star-forming systems with a Icr confidence in- 
terval of 84%-93% and 91%-98%, respectively. The boundaries of the 
four regions are indicated by lines of the same slope but different offset: 
log([NeIII]/[NeII])=8.41og([SIII]/[SiII]+7, where 7=(+3.3,+l. 2,-2.5) for the 
dotted, solid and dashed lines. We include the BCDs (diamonds), starbursts 
(squares) and AGN (crosses) of our sample in the diagram. The coloring 
scheme follows the one of Figures 3 and 5, so that galaxies located at the 
upper and lower branch of Figure 4 are indicated here in orange and cyan 
respectively. Starbursts with a [NeV] 14.32^m detection are indicated with a 
circle around their corresponding symbols. Filled symbols are WR galaxies. 

1999) and/or shock excitation (Lutz et al. 1998). We have al- 
ready shown in §3.2 and §3.3 that the [OIV] diagram can suc- 
cessfully identify galaxies where an AGN contributes to their 
mid-lR spectrum as they are placed on the upper branch of the 
Figure 1 . In this section, we present a preliminary analysis on 
the nature of the [OIV] emission in the starburst galaxies and 
BCDs of the lower branch, focusing on the interpretation of 
models which include WR stars and shocks. More detailed 
study on this subject will be presented in an upcoming paper. 

Some of the sources in our sample are classified as 
Wolf-Rayet galaxies. They show features of WR stars in 
their optical spectra, most commonly the broad HellA4686 
emission which has an IP of ^ 54eV, similar to [OIV]. 
Therefore, identifying correlations between WR features and 
[OIV] emission could, in principle, reveal the contribution 
of those stars. We cross-matched our BCDs with the WR 
galaxy catalog presented in Schaerer et al. (1999) and found 
that seven BCDs (lZwl8, SBS0335-052E, UM461, llZw40, 
NGC1569, Mrkl450, and NGCl 140) as well as two stai'bursts 
(NGC1614 and NGC7714) have reported WR features. 

Unfortunately, this WR galaxy sample is small, and there 
are no consistent measurements of the optical WR features 
from the literature. Thus a direct correlation of the WR 
strength and the [OIV] line emission could not be established. 
In Figure 1 we have marked these WR galaxies (Schaerer et al. 
1999) as filled symbols. It appears that at least the detec- 
tions of WR signatures are not closely associated with high 
[01V]/[S111] ratios, instead, they appear to be more correlated 
with high [Nelll]/[Nell] ratios. 

Another way to probe the contribution of the WR stars to 
the [OIV] emission is by photo-ionization modeling. We per- 
form a tentative exploration of the contribution of those stars 
to the presence of the O^^^ ion in our lower branch galaxies 
with Stai-burst 99 (Leitherer 1999) and CLOUDY (Ferland 
2001). We use Starburst 99 to model the stellar radiation field 
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of an instantaneous burst of star formation aging from 1 to 10 
Myrs. All stars have a metallicity of 1/3 Z© and the follow a 
Salpeter IMF. The output from Starburst 99 is used as input to 
CLOUDY in order to estimate the flux of the mid-lR emission 
lines. 

The models predict that the presence of considerable [OIV] 
emission is directly related to the appearance of the Wolf- 
Rayet phases, becoming strong around 4 Myrs after the onset 
of the burst and fading to undetectable levels 2 Myrs later. On 
the other hand, the [Nelll] line is related to the presence of 
massive main sequence O stars but also becomes very weak 
after 6 Myrs. Both [Nell] and [Sill] Unes are still strong after 
10 Myrs. Consequently the [01V]/[S11I] and [Nelll]/[NelI] 
line ratios peak around 4-5 Myrs after the onset of the burst 
with values of ^^0.6 and 40 respectively. The value of 0.6 for 
the [01V]/[SI1I] line ratio agrees with the highest value ob- 
served in the galaxies of the lower branch in Figure 1, but the 
predicted [NeIII]/[NelI] ratio is about a factor of 4 higher than 
the observed value. Clearly a more systematic study exploring 
thoroughly the parameter space of the models and comparing 
them to larger sample of galaxies with mid-lR measurements 
is required to further investigate this discrepancy. 

We also study the effects of shocks on the production of 
[OIV] utilizing infrared lines which are strong in regions of 
shocks, such as the [Fell]25.99/im (see Lutz et al. 2003). 
The infrared [Fell] lines have been suggested to trace well 
of the supernova activity (Greenhouse et al. 1991, see also 
O'Halloran et al. 2008). The [Fell] ionic abundances in su- 
pernova remnants are orders of magnitudes larger than the 
typical values found in Orion-like Hll regions (e.g. Seward 
et al. 1983; Graham et al. 1987; Oliva et al. 1989). This 
is primarily due to shock destruction of grains. Since iron 
is normally depleted onto the grains, the abundance of iron 
in the ISM is significantly increased in supernova remnants. 
Furthermore, the cooler and denser gas behind a fast shock 
generates more low-ionized Fe^ compared to what is found 
in typical Hll regions. Following the analysis of Lutz et al. 
(2003) to evaluate the effects of shocks in NGC6240, we plot 
in Figure 7 the ratio of [Fell]25.99^m/[01V]25.89^m as a 
function of [01V]25.89/im/[Slll]33.48/im. Also plotted are 
the measurements of the supernovae remnant RCW103 (Oliva 
et al. 1999). Most BCDs and starbursts in our sample show an 
anti-correlation on the diagram, with no obvious indication of 
enhanced [Fell] emission. However, one galaxy, lZwl8, ap- 
pears to be an outlier compared to other BCDs with a clear en- 
hancement of the [Fell] emission relative to the [OIV] emis- 
sion. This enhanced [Fell] emission could indicate a signifi- 
cant contribution of shocks in IZwl8. 

The influence of shocks to the observed [OIV] emission can 
be better constrained with detailed shock models (e.g. Allen 
et al. 2008) but we defer this work to a future paper on this 
topic. Overall, our preliminary study described in this sec- 
tion suggests that there is no strong indication that shocks 
play an important role in the mid-lR emission of starbursts 
and BCDs. Hence, WR stars are the most likely culprit of the 
[OIV] emission seen in the lower-branch galaxies of the Fig- 
ure 1 . We do note though the difficulty to simultaneously fit 
the [01V]/[S111] and [Nelll]/[Nell] ratios with cuiTent photo- 
ionization models. It appears that the models require mecha- 
nisms to significantly suppress [Nelll]/[Nell] ratios for a fixed 
[OIV] emission. This issue is discussed more in §4. 

3.5. The [NeV] 14.32 fi}7i Emission 
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Fig. 7.— A plot of the [OIV]25.89/im/[SIII]33.48/im as a function of the 
[FeII]25.99/xm/[OIV]25.89Atm ratio for the BCDs (diamonds) and starbursts 
(squares). Galaxies located on the upper and lower branch of Figure 4 are 
plotted in orange and cyan respectively. The large asterisk indicated the su- 
pernovae remnant RCW03 (Oliva et al. 1999), while the filled symbols are 
WR galaxies. 

As we have discussed earlier, BCDs display a signif- 
icantly harder radiation field than typical starbursts with 
high [01V]/[S1I1] and [Nelll]/[Nell] ratios. However, no 
[NeV]14.32^m emission has been detected in any of them. 
It is known that [NeV], with an IP of 97eV, serves as an un- 
ambiguous indicator of an AGN and the 7 out of our 24 star- 
burst galaxies where [NeV] 14.32/im was detected, all host an 
AGN (see table 1 of Bernard-Salas et al. 2009). In Figure 8, 
we plot the [NeV]14.32/xm/[Nell]12.81^m ratio as a func- 
tion of [Nelll]15.56/im/[Nell]12.81^mfor BCDs, stai'bursts, 
AGNs, and ULlRGs and find a similar bifurcation signature to 
the one in Figure 1 . This also confirms the bifurcation nature 
of the [OIV] diagram discussed in §3.2 and §3.3. A similar 
bifurcation is also seen in the [NeV]14.32/im/[Nell]12.81/im 
vs. [SlV]10.512Aim/[Slll]18.712/xm diagram (not plotted). 
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Fig. 8.— A plot of the [NeV]14.32^tm/[NeII]12.81/im as a function of 
the [NeIII]15.56^m/[NeII]12.81/tm ratio for BCDs, starbursts, AGNs, and 
ULlRGs. The symbols are the same as in Figure 1. Note that since no [NeV] 
was detected in any BCD only upper limits are indicated for the correspond- 
ing line flux ratios. 



4. GENERAL EXCITATION OF STAR FORMING GALAXIES 

In addition to the high-ionization emission lines 
of [NeV] and [OIV], other common mid-lR emis- 
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sion lines such as [SIV]10.51^m, [SIII]18.71/im, 
[NeII]12.81^m, and [NeIII]15.56^m have been used 
in order to probe the ionization field in BCDs and 
starbursts, the [NeIII]15.56^m/[NeII]12.81/im vs. 
[SIV]10.51^m/[SIII]18.71^m diagram being such an 
example (e.g. Verma et al. 2003; Dale et al. 2006; Farrah et al. 
2007). Both Wu et al. (2006) and Bernard-Salas et al. (2009) 
studied this diagram independently for different samples. 
Here we combine both and also include the data from our 
comparison sample of AGN, ULIRGs, and HII regions. The 
study is similar to the one presented by Groves et al. (2008), 
who also include samples of various galaxy types. The 
advantage of the present work is that the measurements are 
more consistent as all data were obtained with the Spitzer 
high-resolution modules. 
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Fig. 9.— a plot of [NeIII]/[NeII] vs. [SIV]/[SIII] for BCDs, starbursts, 
AGNs, Galactic/LMC/SMC HII regions and ULIRGs. The symbols of BCDs, 
starbursts, AGNs, and ULIRGs are the same as in Figure 1 . The HII regions 
are designated with triangles. Most BCDs, starburst galaxies, and HII regions 
are distributed in a narrow band, indicated with two dashed lines. 



In Figure 9, we plot the [NeIII]15.56Aim/[NeII]12.81/im 
and the [SIV]10.51^m/[SIII]18.71^mratios for the starbursts, 
BCDs, AGNs, HII regions, and low-redshift ULIRGs. For the 
HII regions and the ULIRG sample, only sources with de- 
tections in all four lines are displayed. In addition, we in- 
clude the SH mapping observations of the central positions of 
NGC5253 (Beirao et al. 2006), a typical BCD with metallic- 
ity of 1/6 Zq (Kobulnicky et al. 1999). Since WR features 
have also been found in the center of NGC5253, we desig- 
nate it with a filled diamond as other BCDs with detected WR 
features. 

Observing the plot we note that starbursts populate the 
lower-left corner of the plot (see also Bernard-Salas et al. 
2009). BCDs in general have higher values in both line ra- 
tios than the starbursts, while those with WR signatures sys- 
tematically have even higher [NeIII]/[NeII] and [SIV]/[SIII] 
ratios than those without. The values of [NeIII]/[NeII] in 
BCDs reach as high as 20, and [SIV]/[SIII] as high as 10, val- 
ues which in some cases are even greater than those seen in 
AGNs. BCDs, starbursts as well as giant HII regions all lie in 
a band stretching over 3 orders of magnitude in [NeIII]/[NeII] 
and [SIV]/[SIII]. The band, denoted by the two dashed lines, 
has a slope of ^ 0.9, similar to the slope obtained by Groves 
et al. (2008). As discussed by Groves et al. (2008), AGNs and 
ULIRGs are offset from the band displaying slightly higher 
[SIV]/[SIII] ratios relative to their [NeIII]/[NeII] ratios. This 



is very likely due to the harder radiation field in AGNs which 
excites [SIV] more easily than [Nelll]. However, we notice 
that the offset is not as significant as in the [OIV] diagram. 
This indicates that the [NeIII]/[NeII] and [SIV]/[SIII] ratios 
are not very sensitive to the excitation mechanisms, and can- 
not be used as diagnostics of the central energy source in a 
galaxy (Lebouteiller et al. 2008). 

The higher [NeIII]/[NeII] and [SIV]/[SIII] ratios in BCDs 
compared to average starburst galaxies indicate a harder ra- 
diation field in BCDs. This had already been well estab- 
lished (e.g. Verma et al. 2003; Dale et al. 2006, and ref- 
erences therein). We note, though, that in our sample of 
BCDs and starbursts there is little overlap in the values of 
the [NeIII]/[NeII] and [SIV]/[SIII] ratios with a separatrix at 
[NeIII]/[NeII]- 1.2 and [SIV]/[SIII]- 0.4. Since our BCDs 
have preferentially sub-solar metallicities** (Wu et al. 2008a), 
while the galaxies in the starburst sample are more metal rich 
(Bernard-Salas et al. 2009), their separation in the diagram 
suggests that metallicity is likely the dominant parameter de- 
termining the [NeIII]/[NeII] and [SIV]/[SIII] ratios. 
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Fig. 10. — Top: A plot of the [NeIII]/[NeII] ratio as a function of the Ne 
abundance relative to hydrogen for BCDs, HII regions and starbursts with 
available abundance measurements. Bottom: Similar plot for the [SIV]/[SIII] 
ratio relative to the the S abundance. The object with the lowest measured 
Ne and S abundance is IZwl8 and Tol 65 (Wu et al. 2008a). In both plots the 
vertical dashed line indicates the corresponding solar abundance (see text). 

Here we adopt (Ne/H)Q=1.2x 10"^ from Feldman & Widing (2003), 
(S/H)0=1.4x 10-5 f^^^ Asplund et al. (2005) and (O/H)0=4.9x lO""' from 
Allende Prieto et al. (2001). 
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To examine this more directly, we use the metallicity mea- 
surements of BCDs, HII regions, and starbursts from Wu et al. 
(2008a), Lebouteiller et al. (2008), and Bernard-Salas et si. 
(2009), and plot in Figure 10, the [NeIII]/[NeII] ratios vs. the 
Ne abundances and [SIV]/[SIII] ratios vs. the S abundances. 
We observe an anti-correlation between excitation and abun- 
dances for both Ne and S, even though there is a substantial 
scatter (see also Verma et al. 2003; Wu et al. 2006). 

As mentioned earlier, the low value of the [NeIII]/[NeII] ra- 
tio in starbursts with high metallicities was noted already us- 
ing ISO data (see also Thornley et al. 2000; Verma et al. 2003; 
Madden et al. 2006), but its implications are still puzzling. 
The low ratios apparently suggest a lack of massive stars 
in the starbursts, in contradiction to the fact that high-mass 
stars are clearly observed in nearby starbursts (e.g. O'Connell 
et al. 1995; Whitmore & Schweizer 1995; Meurer et al. 1995). 
Thornley et al. (2000) attributed this to an aging effect: mas- 
sive starbursts have a shorter lifetime, thus, the hottest stars 
may not be the dominant source of ionization in the star- 
bursts. However, Rigby & Rieke (2004) propose that as high- 
mass stars spend much of their lives embedded within ultra- 
compact HII regions (i.e. Charmandaris et al. 2008), the high 
density and high extinction of the ultra-compact HII regions 
prevents the mid-infrared nebular lines from forming and es- 
caping. 

We note that the starbursts with low [NeIII]/[NeIIJ ratios 
have fairly high [OIV]/[SIII] ratios (see Figure 1). This puts 
additional constraints to the issue of low values of observed 
[Neni]/[NeII] ratios. As discussed in §3.4, if WR stars are 
the main source of the [OIV] emission, the observed values 
of the [OIV]/[SIII] in starbursts argue strongly against a low 
mass cutoff of the IMF in these starbursts, since the WR phase 
would be too short to sustain the intensity of the [OIV] ob- 
served. It will also places some constraints on the scenario 
that massive stars spend much of their life in ultra-compact 
HII regions. The strong extinction and high density will have 
to somehow significantly affect [Nelll] and [Nell], but not as 
much the [OIV] emission. 

In Figure 10, we notice that the two extremely low metal- 
Ucity objects, IZwlS and Tol65 (see Wu et al. 2008b), de- 
viate from the anti-correlation of the emission line ratios vs. 
the abundances. It appears as if the anti-correlation starts af- 
ter a maximum value for [NeIII]/[NeII] and [SIV]/[SIII] at 
~20 and ~8 respectively. We should note though that ex- 
tremely low metallicity galaxies are scarce, and their abun- 
dance measurements using the infrared fine-structure lines 
are more challenging due to the non-detection of their Hua 
I2.7/Ltm fine. However, similar plots using the oxygen abun- 
dance determined from optical spectroscopy also shows the 
same scatter (Wu et al. 2006). 

5. CONCLUSIONS 

We have presented our analysis of the high spectral resolu- 
tion mid-IR observations for a sample of 12 BCDs obtained 



with Spitzer/IRS. A direct comparison with samples of star- 
burst galaxies, AGN, ULIRGs, as well as isolated HII regions, 
enabled us to extend earlier studies of the ionization and ex- 
citation in star forming galaxies towards lower metallicities, 
and to use the observed mid-IR line emission to explore the 
excitation mechanisms in these systems. 

Although BCDs are known to have a hard radiation field, no 
[NeV]14.32/im emission has been detected with Spitzer/IRS 
even for the most metal-poor BCDs. 

We find that the [OIV]25.89/im line is commonly detected 
in the BCDs, starbursts and AGNs. We have demonstrated 
that the diagram of fOIV125.89/im/[SIII]33.48/im as a func- 
tion of |NeIII|15.56^m/[NeIIJ12.81/im can be useful diag- 
nostic to probe the dominant power source in dust enshrouded 
systems. Star forming galaxies populate two branches in 
this diagram, with those harboring an AGN component lo- 
cated on the upper branch, and those without on the lower 
one. This diagram provides a classification consistent with 
the log([NII]/HQ:) vs. log([OIII]/H/3) diagram in the opti- 
cal, especially for [NeIII]/[NeII] values above 0.3. We em- 
pirically define a line (Equation 1) separating galaxies with 
and without AGN contribution, which could be used in con- 
junction with other diagnostics to better determine the pres- 
ence of an AGN in the infrared. A similar bifurcation is 
seen in a [NeV]/[NeII] vs. [NeIII]/[NeII] diagram. In the 
[NeIII]/[NeII] vs. [SIV]/[SIII] diagram. AGNs and ULIRGs 
are only slightly off the band defined by BCDs, starbursts, 
and HII regions, with higher [SIV]/[SIII] values for given 
[NeIII]/[NeII] ratios. 

In an effort to probe the contribution of shock excitation 
to the observed [OIV]25.89/xm flux, we compare it with the 
[FeII]25.99/im line , which is a good tracer of supernovae 
activity. We found no evidence of shock contribution, demon- 
strated by an enhanced [Fell] emission, in starbursts or BCDs. 
The only exception is IZw 18, which may indicate a significant 
shock contribution in this low metallicity system. 

The most likely contributor to the [OIV] flux, both in star- 
bursts and BCDs, are WR stars. However, current photo- 
ionization models predict flux ratios of [NeIII]/[NeII] which 
are too high for the observed [OIV]/[SIII] values. This places 
an additional constrain on the problem of low [NeIII]/[NeII] 
ratios already known in high-metallicity starbursts, and will 
be further investigated in future work. 
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Fig. 1. — Zoom in of IRS high resolution spectra near the wavelength range of [SIV]10.51/im, [NeII]12.81/im, [NeIII]15.56/im, and [SIII]18.71/im for 
NGCl 140, NGC1569, IIZw40, IZwl8, Mrkl499, and Harol 1. The expected location of the each emission Une is indicated with the dashed vertical line. 
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Fig. 2. — Zoom in of IRS high resolution spectra near the wavelength range of [SIV]10.51^tra, [Nell] 12.8 l/.(m, [Nelll] 15.56^m, and [SIII] 18.71/im for 
SBS0335-052E, UGC4274, Mrkl450, UM461, Tol65, and VIIZw403. The expected location of the each emission line is indicated with the dashed vertical line. 



13 



0.066 ■ 
0.064 ; 
0.062 - 
0.060 ; 
0.058 ■ 
0.056 i 



[NeV]14.3^m 

NGCl 1 40 



[0IV]25.9^m 

NGCl 14i 



[Slll]33.5/xm [Sill]34.8/xm 




. NGCl 14| 

A 


1 ' ^ 


. NCC156'] 




" IIZw40 1 




IZw18 


' ' 1 ' i ' 1 i ' 


Mrk1499 


' ' 1 ' i ' 1 i ' ' 


- Harol 1 j 





NGCl 1 4o[ 


1 


NGG1569 
^""''i''^V7""r~!i 


' ' 1 ' ' i 1. 


- IIZw40 1 


' ' ' ' 1 ' ' i 


IZw18 / 
1 i ' ' 1 ' ' i 1' 


' ' ' 1 ' ' i 1 
' ' ' 1 ' " i 


Mrk1499 




" Harol 1 1 





14.0 14.2 14.4 14.6 14.8 



25.6 25.8 26.0 26.2 26.4 33.0 33.2 33.4 33.6 33.E 



34.4 34.6 34.8 35.0 35.2 



Rest Wavelength(/xm 



Fig. 3. — Zoom in of IRS high resolution spectra near the wavelength range of [NeV]14.32^m ( with its adjacent [CI II] 14.37^tm), [OIV]25.89^m, 
[FeII]25.99/im, [Sin]33.48/im, and [SiII]34.8I/im forNGCIMO, NGC1569, IIZw40, IZwlS, MrkI499, andHarolI. The expected location of the each emission 
Une is indicated with the dashed vertical Une. 
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Fig. 4. — Zoom in of IRS high resolution spectra near the wavelength range of [NeV]14.32/im (with its adjacent [CI II]14.37^m), [OIV]25.89^m, 
[FeII]25.99/im, [Sin]33.48/im, and [SiII]34.81/im for SBS0335-052E, UGC4274, Mrkl450, UM461, Toll214-277, Tol65, and VIIZw403. The expected 
location of the each emission line is indicated with the dashed vertical line. 



TABLE 1 

Mid-Infrared Line Measurements of BCDs 



HuxMxlO-i^Wni-^) 

Object Name [SIV] [NeH] [NelH] [Sm] [SIII] [Sill] [NeV] [NeV] [OIV] [Fell] 

10.51//m 12.81/im 15.56;um 18.71/im 33.48/im 34.81/im 14.32/im 24.32;im 25.89/im 25.99;im 

34.8eV 21.6eV 4L0eV 23.3eV 23.3eV 8.2eV 97.1eV 97.1eV 54.9eV 7.9eV 
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"The integrated fluxes are measured iirom the liigh resolution spectra of the targets. Background emission has not been subtracted and no scaling between the SH and LH modules of IRS has been applied. 
The central wavelength and the ionization potential of each line are indicated. 



